Singlet Ground State of the Quantum Antiferromagnet Ba 3 CuSb 2 09 
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We present local probe results on the honeycomb lattice antiferromagnet BaaCuSb2 09. Muon 
spin relaxation measurements in zero field down to 20 mK show unequivocally that there is a total 
absence of spin freezing in the ground state. Sb NMR measurements allow us to track the intrinsic 
susceptibility of the lattice, which shows a maximum at around 55 K and drops to zero in the low- 
temperature limit. The spin-lattice relaxation rate shows two characteristic energy scales, including 
a field-dependent crossover to exponential low-temperature behavior, implying gapped magnetic 
excitations. 



In recent years, the field of geometrically frustrated 
magnetism has been galvanized by the discovery of sev- 
eral promising candidates for quantum spin liquid (QSL) 
physics in kagome [TJ [2] , hyperkagome [3] and triangular- 
lattice [H [5] geometries. The most recent and surprising 
discoveries in this vein have come in the 6H-perovskite 
family, Ba 3 MSb 2 09 where M =Cu, Ni [3 [7], and have 
generated an enormous theoretical interest [8Ul2j. 

Here we focus on the S = 1/2 antiferromagnet 
Ba3CuSb2 0g, originally thought to feature a triangular 
lattice of Cu 2+ spins [13] and shown by Zhou et al. to 
exhibit a very surprising gapless QSL state [B] with a 
lack of static magnetism down to 200 mK, well below the 
Weiss temperature, 9w — 50 K, and a peculiar linear spe- 
cific heat. An in-depth characterization, including single 
crystals, of Nakatsuji et al. 14] on the other hand has 
shown that the structure more likely involves triangular 
bilayers equally occupied by Cu and Sb that form a deco- 
rated honeycomb lattice during crystal growth as a result 
of electric-dipole interactions between Sb 5+ -Cu 2+ 'dumb- 
bells' that map onto a frustrated Ising model. Moreover, 
it has been shown that a crucial ingredient in the physics 
of Ba 3 CuSb 2 09 is the orbital degrees of freedom. The 
system is Jahn- Teller (JT) active, with three equivalent 
local distortions and corresponding d x 2_ y 2 orbitals. In 
some off-stoichiometric samples, a collective orthorhom- 
bic JT transition at ~ 200 K occurs, whereas stoichiomet- 
ric samples maintain hexagonal symmetry, show minimal 
anisotropy and remain dynamic. 

The lack of static magnetism in Ba3CuSb20g to T <C 
®w ~ 50 K remains a profound result. Neither trian- 
gular nor honeycomb nearest neighbor antiferromagnets 
are expected to have QSL ground states and apparent 
JT instabilities make it even more surprising. Two likely 
proposals have emerged: either a new type of QSL on 
the honeycomb lattice accompanied by dynamic JT dis- 
tortions or a random singlet state (RSS) driven by static 



JT distortions on a rather disordered lattice O |T5] . An 
understanding of this system therefore hinges on a char- 
acterization of the ground state and the evolution of the 
susceptibility and spin fluctuations. 

In this Letter, we have applied the local probe tech- 
niques, muon spin relaxation (/xSR) and nuclear magnetic 
resonance (NMR), to Ba3CuSb2 0<). We demonstrate in 
this system a lack of static magnetism in zero-field, sin- 
glet formation around 50 K and a lower temperature 
field-induced gap. Our measurements help distinguish 
between the two possible scenarios proposed in Ref. [14] 
and favor a RSS. 

The synthesis and structural characterization of the 
polycrystalline Ba3CuSb2 0g samples studied in this work 
are described in the Supplemental Material (SM) [16] , 
Macroscopic susceptibility, Xmacro, measurements reveal 
a Weiss temperature, 9\y ~ 51 K, consistent with previ- 
ous work [Bldl] and Curie constant 0.48 Kemu/mol/Oe, 
implying either a somewhat large g = 2.27(1) (relative to 
ESR results, g = 2.20 Q3]) or else a 6% surplus of Cu. All 
characterization techniques employed [16] have revealed 
an Sb-Cu ratio stoichiometric to well within ±10% indi- 
cating [T3] that our samples should maintain hexagonal 
symmetry. 

To further confirm that our samples are consistent with 
the hexagonal phase of Ref. |il4] and not the orthorhom- 
bic phase that exhibits freezing at T g = 100 mK, we 
performed fiSK experiments at PSI from 150 K down to 
20 mK. The zero field (ZF) polarization of the bulk [17] 
can be fit at all temperatures with a static Kubo-Toyabe 
(KT) function with AH = 0.76 ± 0.2 G, resulting from 
nuclear moments, and weak exponential relaxation of dy- 
namical origin (see Fig. [I]). The results show that there 
is a complete lack of spin freezing [TH] in BaaCuSb20g 
down to 20 mK in ZF, in agreement with Ref. [14]. Ap- 
plying a 50 G longitudinal field (LF) decouples the muons 
from the nuclear moments and reveals changes in spin 
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FIG. 1: (a) Muon polarization at T — 10 K in zero field and 
at T = 20 mK in zero field, and longitudinal field of 50 G. The 
fits are exponential relaxation, multiplied with a KT form of 
nuclear origin for the ZF data, (b) Relaxation rate A(T) in 
a longitudinal field of 50 G. The shaded region is a guide to 
the eye 



dynamics. The relaxation rate, A(T), obtained through 
exponential fits (inset of Fig. [I]), is very weak with a sub- 
tle increase as T is decreased and is comparable to that 
of herbertsmithite [TJ [TD] . 

To precisely study the spin dynamics and in-field be- 
havior of Ba3CuSb2 0g, we turn to a more strongly- 
coupled probe: Sb NMR. Spectra are obtained with the 
standard tt/2-t-it pulse sequence ( 14 < r < 23 fj,s) and 
integration of the echo as the field is swept. Below 20 
K, Cu NMR signals are seen, superimposed on the 121 Sb 
spectra. A 135 Ba quadrupolar satellite is found near the 
123 Sb line, but is small and does not change in temper- 
ature. Hence, most of this work concentrates on 123 Sb. 
The high-tcmpcrature NMR spectra are consistent with a 
distribution of quadrupolar couplings, since the broaden- 
ing of the quadrupolar satellites is larger than that of the 
central line. This can be seen to result from Cu-Sb dumb- 
bell disorder inherent in the decorated honeycomb lattice 
and random JT distortions [14.. A powder simulation 
with parameters v Q ~ 3.7±0.9 MHz for the I = 7/2 123 Sb 
nucleus is shown in Fig. [2] While there are two inequiva- 
lent Sb sites in the structure, they are not discernible in 
the spectra. The lineshift, K, obtained from the position 
of the maximum of each spectrum and is shown in Fig. [3] 
A scaling of K at high T with Xmacro reveals a hypcrfinc 
coupling of A = dK/dXmacvo — 12 kOe/^i B and shows 
that the temperature-independent shift (which may be 
chemical and/or quadrupolar) is only ~ 80 ppm. 

The most striking observation in the Sb NMR spec- 
tra is the non-monotonic shift (intrinsic susceptibility), 
K (xint)i shown in Fig. [3} The maximum in Xint is found 
at a remarkably high temperature ~ 55 K ~ 9w- Clearly 
this material does not behave like a nearest-neighbor- 
s' = 1/2 triangular lattice antiferromagnet, as seen from 
comparison with the high-temperature series expansion 
(HTSE) [2D], in Fig. |] The appearance of strong spin 
correlations or dimerization at such a high temperature 
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FIG. 2: (a) Selected 123 Sb NMR spectra, including a sim- 
ulation of the spectrum at 80 K. The 135 Ba satellite can be 
isolated (shaded area) by roughly doubling the rf pulses, (b) 
Peak of spectra at several temperatures showing the non- 
monotonic shift (green circles). 



is unusual for a QSL and implies a lack of magnetic frus- 
tration. 

Equally important is the observation that Xint drops 
out appreciably and appears to approach zero as T — » 0. 
We cannot, however, differentiate between an exponen- 
tial (gapped excitations) or power-law (gapless excita- 
tions) limit with K given the sizable linewidth. Xint here 
contrasts heavily with behavior observed in other QSL 
candidates. For example the organic triangular lattice 
systems, which are thought to be described by a QSL 
model with a Fermi surface of spinonspT] [22] . exhibit 
finite susceptibility as T — > [HE]. 

Comparison with Xmacro (Fig. [3]) illustrates the signif- 
icant level of defect contribution in the system. Magne- 
tization measurements at low T up to 14 T, show that 
~ 16% of the spins in the system are weakly interacting 
with #dcf — —0.5 K at the lowest temperatures measured 
and are easily saturated with magnetic field [15] . Note 
that the sample studied here displays the same size of 
Curie tail as in Refs. [DJ [H] suggesting that it might be 
an intrinsic property of the decorated honeycomb lat- 
tice. Although such a weakly coupled defect or 'orphan 
spin' [24 contribution is commonplace in related sys- 
tems [23] [251427] , here it seems likely that it originates 
from the out-of-plane Cu' site which is inherent in the 
decorated honeycomb structure and is linked to the hon- 
eycomb motifs by a frustrated isosceles triangle (see the 
structure in Ref. [T4"]L 

Meanwhile, the linewidth (shown in the inset of Fig. [3j 
obtained with Gaussian fits to parts of the spectra where 
the central line can be isolated) monotonically increases 
despite the drop in Xint, implying that the linewidth is 
dominated by a defect-induced staggered magnetization 
or spin texture related to the Curie-tail. At low T, the 
linewidth surprisingly represents a significant fraction of 
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FIG. 3: Left axis: NMR lineshift, K, taken with 121 Sb at 
6.4 T (high T only) and 123 Sb at 12 T. Right axis: Xmacro 
at 5 T minus diamagnetic susceptibility Xdia = —1 x 10~ 4 
emu/mol/Oe [H]. Also shown is the HTSE of \ for the S = 
1/2 triangular lattice antiferromagnet [20]. Inset: Standard 
deviation of local magnetization, a m = AAH. 
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FIG. 4: (a) 123 Sb spin-lattice relaxation rate, l/Ti(T), taken 
at 6.5 T and 12 T. (b) Semilog plot of 1/Ti vs. 1/T showing 
gapped behavior at low temperature, with A ~ 6 K at 6.5 T 
and A ~ 10 K at 12 T. Inset: the recovery curves M(t) at 
selected temperatures showing a broadening distribution of 
relaxation times at low T, that is nonetheless small relative 
to the orders of magnitude drop in 1/Ti. 



the full moment (almost 0.4 hb), despite zero lineshift. 

To assess the nature of the magnetic excitations, we 
turn to the spin-lattice relaxation rate, 1/Ti. Because 
the recovery curves, M(t), consist of multiple stretched 
exponentials [28] , we have simply defined T± by the point 
where 1 — M(t)/Moo = 1/e. The quadrupolar Sb nuclei 
are in principle sensitive to magnetic and charge fluctua- 
tions. Comparison of 1/Ti for the two Sb isotopes (with 
very different gyromagnetic ratios, 7, and quadrupolar 
moments, Q) at constant jH and T > 20 K, shows a ra- 
tio 123 Ti/ 121 Ti ~ 1.6. This is close to the ratio 1.7 that 
we would expect if we were primarily probing magnetic 
fluctuations, rather than 0.6 which would be expected of 
charge fluctuations [28 . Thus the JT distortions pro- 
posed in Ref. [Uj are either fully static or dynamic on a 
time scale too fast for NMR. 

The T-dependence of 1/Ti, taken at fields of 6.5 T and 
12 T is shown in Fig. [4] and shows three distinct regimes 
of relaxation. At high T, a small drop in 1/Ti is ob- 
served and in a plot of log(l/Ti) vs. 1/T (Fig. gb)) 
can be linked to a ~ 50 K gap. This, along with the 
suppression of x starting ~ 55 K, seems to show singlet 
formation, with an energy scale of Ai ~ J ~ 9w — 50 K 
and is consistent with a broad peak in the inelastic neu- 
tron scattering (INS) spectra of Nakatsuji et al, which 
was attributed to the formation of short-range singlets on 
hexagonal motifs [14]. However, that drop in relaxation is 
quickly overtaken by appreciable low energy fluctuations 
giving a largely temperature independent relaxation from 
5 to 20 K. This quasi-plateau in 1/Ti is likely related to 
the continuum seen in INS [T3] . 

There is another, important energy scale that is driven 



by the applied field. Below the quasi-plateau, a crossover 
to an exponential drop occurs, indicating a true gap to 
magnetic excitations rather than the partial dimerization 
at ~ 50 K. The crossover temperature, Tc, is found to 
be roughly proportional to the applied field: Tc — 5 K 
at 12 T and T c ~ 2.6 K at 6.5 T. The inset of Fig. g) 
shows the expected T-dependence, 1/Ti oc exp(— A /T), 
of gapped excitations and a field-dependent gap A (H). 
At 12 T, A ~ 10 K and at 6.5 T, A ~ 6 K. The Ti 
recovery curves (inset of Fig.|4|, show that below ~ 10 K, 
the distribution of relaxation times widens appreciably, 
suggesting that there is in fact a distribution of gaps, as 
might be seen in a RSS 29 . 

This behavior runs contrary to the case of a standard 
gapped dimer state or a gapped QSL where a magnetic 
field closes the energy gap and eventually leads to mag- 
netic order and power-law behavior [301 IS]- It is also 
surprising given the C(T) of Zhou et al. [61 which (after 
subtraction Schottky anomalies) is very small but linear 
in the low-T limit and field-independent to as high as 9 
T. This is difficult to reconcile with our observation of a 
clear gap to magnetic excitations, although a theoretical 
model supporting such a dichotomy between C oc T and 
exponential 1 jT\ has been proposed for the related 5=1 
system Ba 3 NiSb 2 9 [8]. 

On the basis of the structural and orbital informa- 
tion obtained by Nakatsuji et al., two main scenarios for 
Ba 3 CuSb 2 9 have been suggested [H] [T5]: (1) A dy- 
namical JT distortion that allows the system to main- 
tain isotropic character may lead to an exotic spin-orbital 
liquid state with gapless excitations. (2) A disordered 
static JT distortion may occur that stabilizes a RSS. The 
configuration of distortions was proposed to obey locally 
a three-fold rotation axis, with three dimers around a 
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hexagon of Cu atoms. 

If the first scenario is correct, the gap seen in 1/Ti 
likely represents a transition to a frozen state with 
gapped magnons. Indeed our 1/Ti results show qualita- 
tive similarities to the relaxation in other gapless QSL 
candidates [3"TH3~4"] which exhibit field-induced phases 
with steep power law or exponential l/Ti(T). In all cases, 
the change in dynamics is a smooth crossover, with no 
peak in 1/Ti, atypical of standard phase transitions. A 
freezing of spins should be accompanied by a change in 
linewidth, but in this instance, it is difficult to distinguish 
between defect induced broadening and signs of freezing. 
Much of the increase in linewidth (Fig. [3]) is occurring 
well above the crossover temperature (~ 5 K at 12 T) 
and is likely related to the defect contribution in Xmacro- 

However, the QSL scenario (1) is put in doubt by 
strong evidence of dimerization occurring at ~ 55 K, 
as high as the exchange energy J. This is inconsistent 
with QSL phenomenology, where gaps are either absent 
or much smaller than J [35j . and more befitting of an 
unfrustrated dimer system (35] ■ Given our T\ and Xint 
results, a more likely scenario is a type of RSS (2) with 
energy gap of Ai ~ 50 K. As proposed by Nakatsuji et 
al. [2], this RSS consists largely of short-range singlets 
and our work is not consistent with a state involving ar- 
bitrarily long-range singlets ^37 that would result in a 
finite susceptibility at low T. 

Clearly, however there remains a continuum of ex- 
citations below Ai and the field- induced gap A is 
much smaller. Similarity of the field-dependent gap size, 
Aq/H ~ 0.8 K/T, with the Zeeman energy for an isolated 
5 = 1/2 moment, E z /H ~ 1.34 K/T, leads us to pro- 
pose that it may be the weakly interacting defect spins 
that give rise to low-energy excitations in zero-field and 
that freeze out as the defects are saturated by magnetic 
field. Specific heat measurements are consistent with this 
picture, showing Schottky anomalies that shift to higher 
T with applied field [ICS]. 

The effect of these defect spins on the Sb nuclei is likely 
indirect since no shift is seen at low T and the relaxation 
in the plateau regime is too fast to be coming from dis- 
tant spins. With a weak coupling between defects, so 
6*d f — —0.5 K, these defects would be slowly fluctuating, 
but even so, if they are coupled to the nuclei with the 
dipolar interaction, 1/Ti would be an order of magnitude 
smaller than what is seen. It is therefore more likely that 
they are distributed defects that cause a staggered mag- 
netization or spin texture which is then strongly coupled 
to the Sb nuclei. As such, the spin textures contribute 
heavily to the NMR linewidth, but negligibly to K. This 
is commonplace in frustrated [231 [27] and gapped quan- 
tum spin systems [38] . The precise nature of defects and 
resulting spin textures remains difficult to assess but they 
likely result from the out-of-plane Cu' site which is linked 
to the honeycomb motifs by a frustrated isosceles trian- 
gle (see the structure in Ref . [2] ) ■ The Cu' spins should 



be strongly coupled to the lattice [13], but may become 
unconstrained and easily saturated under applied field as 
the lattice becomes dimerized at low temperatures. 

To conclude, we have performed local probe measure- 
ments on the 5 = 1/2 honeycomb system Ba 3 CuSb 2 09, 
revealing an exotic non-magnetic state in zero-field with 
dimerization occurring at 50 K and a lower temperature 
field-induced gap to magnetic excitations. Our measure- 
ments primarily support a RSS with low-energy excita- 
tions resulting from defect-induced spin textures that can 
be saturated with magnetic field. This work allows us to 
distinguish between two proposed [TH [T3] orbital scenar- 
ios and likely implies random but static JT distortions. 
Within this new context of short range honeycomb struc- 
tural order and orbital degrees of freedom, it would be 
enlightening to similarly investigate the 5=1 equiva- 
lent Ba3NiSb20g which also may show a dynamic ground 
state [7J but may not be JT active. 
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SUPPLEMENTAL MATERIAL 



Here we provide additional information on 1) the synthesis and chemical analyses of powder samples of Ba 3 CuSb 2 09, 
2) the analyses of powder X-ray diffraction data, 3) the details of spin-lattice relaxation and 4) high-held susceptibility 
measurements. 



SM.l - Synthesis and chemical analyses of powder samples of BaaCuSt^Og 



Polycrystalline samples of Ba 3 CuSb 2 9 were prepared using a standard solid-state reaction method, starting from 
high-purity BaCC>3, CuO, and Sb2C>3 powder materials. The dried starting materials were examined by x-ray diffrac- 
tion (XRD) measurements. Stoichiometric amounts of these materials were ball-milled, pelletized, and then heated 
at 1100°C for several days in air with several intermediate grindings. Samples were furnace cooled at the end of the 
final heat treatment. As-prepared powders are yellow-brown in color. The samples were analyzed with a scanning 
electron microscope (SEM) equipped with an energy-dispersive spectroscopy (EDS) detector. SEM analyses done in 
both backscattered and secondary electron imaging modes showed the samples to be homogeneous. Semiquantitative 
microprobe elemental EDS analysis was performed at different positions on the sample surfaces using polished sec- 
tions. Atomic ratios of cations were found to be 1.48(3), 2.9(2), and 0.51(3) for Ba/Sb, Ba/Cu, and Cu/Sb ratios, 
respectively. Cation stoichiometries were also determined by means of inductively coupled plasma optical emission 
spectroscopy (ICP-OES) which yielded atomic ratios of 1.44(1), 3.02(3), and 0.48(1) for Ba/Sb, Ba/Cu, and Cu/Sb 
ratios, respectively. These ratios are close to the nominal ones. As the deviation from 2:1 Sb/Cu stoichiometry is well 
below 10%, the results of Nakatsuji et al. [3] suggest that most of the sample should remain in its hexagonal form 
(not undergoing an orthorhombic Jahn- Teller distortion) down to low temperatures. A thermogravimetric study was 
performed in air and in oxygen atmosphere, and no weight decrease or increase was observed below 1200°C. 
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FIG. SM.l: Observed and calculated (solid curve) powder XRD pattern of a sample of Ba3CuSb209 collected at room tem- 
perature using Cu Kai radiation. The tickmarks indicate the positions of the Bragg reflections. The lower curve shows the 
difference between the observed and calculated data on the same scale. 
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SM.2 - Powder X-ray diffraction data for Ba 3 CuSb 2 9 

Laboratory powder XRD data were collected at room temperature in the 29 range of 8-100° using Cu Kai radiation 
(A = 1.5406 A). All patterns showed narrow diffraction peaks and no sign of unreacted starting materials (Note that 
the most intense peaks of Tenorite do not coincide with those of the desired phase). As noticed earlier for other 
Ba 3 MSb 2 0g (M = Mn, Co, Ni) compounds Q], a very small amount of impurity phase BaSb 2 6 (JCDS #82-0520) 
showed up in the diffraction patterns. The data were analyzed with the Rietveld method by use of the JANA 
program [2]. Three 29 regions where the BaSb20g peaks showed up were excluded from refinements. We used 
the recently published PQ^/mmc hexagonal crystal structure [3] as a starting model. The lattice parameters, a = 
5.80390(4) A and c = 14.3231(1) A were in good agreement with those determined from the single-crystal study in 
Ref. [3J. There was no sign of additional structural periodicities. For the Rietveld analysis, independent anisotropic 
atomic displacement parameters (ADPs) for each Ba site, and isotropic ADPs for the other sites, were refined with 
full occupancies of all sites. The two central sites of the face sharing octahedra (Wykoff 4/ positions) were equally 
occupied by Cu and Sb atoms. With this model, patterns were refined with reliability factors of R p = 0.0855, 
Rwp = 0.1225 and GOF = 1.2. Maximum and minimum electron densities in the difference Fourier map were 0.8 
electron/A 3 and -1.2 electron/A 3 , respectively. Empirical bond-valence sums for Bal, Ba2, Sbl, 01, and 02 atoms 
were consistent with nominal oxidation states Ba 2+ , Sb 5+ , 2 ~. We were not able to resolve afew percent change in 
Cu/Sb2 or Cu/Sbl ratios nor a weak oxygen off-stoichiometry. A representative Rietveld profile fit is shown in Fig. 
SMI. 



SM.3 — Isotopic dependence of spin-lattice relaxation 

In Fig. SM2, l/Ti(T) is shown for the two different Sb isotopes at constant frequency jH = 66 MHz. This 
comparison allows us to roughly determine whether our relaxation is dominated by magnetic fluctuations or by charge 
fluctuations, both of which can affect quadrupolar nuclei. The multi-exponent relaxation of a quadrupolar nucleus 
yields different results depending on the spin number I. Here I = 5/2 for the 121 Sb isotope and / = 7/2 for the 123 Sb 
isotope. The parameters of Ref. [I] illustrate that 121 T le / 123 T le ~ 2W\2\/W\2z where W n is the spectral density of 
fluctuating fields for isotope n. Following the calculations of Moriya [S], W n cx 7 2 . Here 7121/7123 = 1-85. Thus if 
fluctuations are magnetic in origin, we expect roughly a factor of 1.7 decrease in 1/T± e as we move from 121 Sb to 
123 Sb. Indeed, this is very close to what is observed at higher temperatures (away from any field induced physics), as 
can be seen from Fig. SM2(a). 

On the other hand, in the case of purely quadrupolar relaxation, we should expect the isotopic dependence to vary 
in the opposite manner as 121 Wq/ 123 Wq ~ ( 123 Q/ 121 Q) 2 = 1.65. This is clearly not the case thus we are primarily 
sensitive to magnetic fluctuations and not to charge fluctuations. Either the charge fluctuations are completely static 
(frozen) in the temperature range of our measurements, or else in the case of a dynamic Jahn- Teller effect [3J, the 
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FIG. SM.2: Ti vs. temperature for the two different isotopes of Sb taken at the same frequency, yH — 66 MHz. 





charge fluctuations may occur at too high a frequency to be probed with NMR. 



SM.4 - High-field magnetization measurements 



High-field magnetization measurements were performed with a Cryogenic vibrating sample magnetometer (VSM) 
system, at temperatures from 2 K up to 5 K. The resulting curves, shown in Fig. SM3, are fit with the formula 

M(H) = n imp (3/2)^ B tanh 

where 9- lmp , and x res are left as free parameters and g is taken from ESR measurements of Ref. [3]. The total number 
of Cu spins in the system, ntot, is taken from a high temperature Curie- Weiss fit using the same <?-factor. The result 
is that n imp /n tot = 0.16, showing that 16% of the spins in the material are easily saturated, with 9 lmp varying from 
-2.4 K (at T = 5 K) to -0.5 K (at T = 2 K). There remains some residual susceptibility, x rcs that is mostly constant 
in T and averages to ~ 2.8 x 10 -3 emu/mol/Oe coming from spins that are not easily saturated. 

Furthermore, we note that there are no features in the magnetization indicative of a phase transition occurring at 
the same field/temperature as the field-induced gap seen in NMR. 



(SM.l) 
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